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We study the dynamical decoupling of multiqubit states from environment. For a system of m 
qubits, the nested Uhrig dynamical decoupling (NUDD) sequence can efficiently suppress generic 
decoherence induced by system-environment interaction to order A*" using (iV + 1) pulses. We 
prove that the NUDD sequence is universal, i.e., it can restore the coherence of m-qubit quantum 
system independent of the details of system-environment interaction. We also construct a general 
mapping between dynamical decoupling problems and discrete quantum walks in certain functional 
spaces. 
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Dynamical decoupling (DD) is a powerful tool to pro- 
tect quantum systems from decoherence induced by the 
inevitable system-environment interaction The idea 
of DD is to dynamically control the system (or environ- 
ment) evolution to suppress the decoherence caused by 
interaction. For example, a static magnetic field with un- 
known magnitude Bz<Jz can induce dephasing of a qubit, 
but such dephasing can be fully eliminated by a spin flip 
ax (i.e., Hahn echo) at half way of the evolution 2J. In 
practice, however, the Hahn echo only suppresses the de- 
phasing to O (T^) for total evolution time T, because the 
magnetic field may have complicated time-dependence in 
both magnitude and orientation due to the evolution of 
the environment. Furthermore, if the environment con- 
sists of quantum degrees of freedom, it can become en- 
tangled with the system via the interaction. Hence it is 
a challenging task to design a universal DD scheme that 
can suppress decoherence to desired order independent 
of the details of system-environment interaction. 

One particular interesting DD scheme is the concate- 
nated DD (ODD), which has been shown to be universal 
for single qubits [3]. The limitation, however, is that 
the pulse number increases exponentially with the sup- 
pression order N (approximately 4^ pulses to suppress 
both bit-flip and dephasing processes to O (T^+^)). It 
is the discovery of the universality of Uhrig DD (UDD) 
sequence [4]-[8] that makes the universal DD practically 
feasible. In contrast to CDD demanding exponentially 
many pulses [3], UDD uses only O (N) spin-flip pulses 
to suppress the dephasing processes to O (T^+^) 
The discovery of UDD sequence has inspired many ex- 
perimental efforts to further improve the coherence over 
a wide range of quantum systems, including trapped ions 
[5] , electron spins [IDj , defect centers [TUl E] , quantum 
dots [Ellin], and superconducting qubits [H]. However, 
UDD is restricted to pure dephasing errors of a single 
qubit. It is desirable to have an efficient DD scheme (with 
poly (N) pulses) to suppress both bit-flip and dephasing 
processes for multiple qubits to O (T^+^). 

Recently, the quadratic DD (QDD) scheme has been 
proposed [15], which uses (A^ + 1)^ pulses to suppress 
both bit-flip and dephasing errors of single qubits. As a 
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FIG. 1: (color online). Nested Uhrig dynamical decoupling 
(NUDD) scheme with 2m nesting levels and suppression order 
N = 2. (a) The timings of NUDD pul ses have a self-similar 
structure, determined by Eqs. ([2| and ( 10 1 . The set of pulses 



associated with rth level is cri for r — 2j — l, and al for r = 
2j. (b) The time dependent modulation functions Fa (t) of 
the toggling frame Hamiltonian, and the corresponding pulses 
for the mth qubit. 



generalization of QDD from 1-qubit system to m-qubit 
system, the nested UDD (NUDD) scheme has been pro- 
posed [T6til8j . which uses (A^ -I- 1)^™ pulses to suppress 
decoherence from the most general interaction between 
the m-qubit system and environment. Although there 
are numerical evidences 15] and theoretical implications 
[711111111] that QDD/NUDD might be universal, it is stiU 
an open question whether QDD/NUDD are universal or 
not [Hl^. 

In this Letter, we shall present a rigorous proof that 
the NUDD scheme with 2m nesting levels and {N + l)^™ 
pulses is a universal DD scheme for m-qubit system, 
which suppresses decoherence processes to O (T^'^^) for 
arbitrary system-environment interaction. We achieve 
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this by providing a mapping between NUDD and a dis- 
crete "quantum walk" in 2m dimensional space. The 
rules that govern this this quantum walk do not depend 
on 2m, which allows for an efficient proof for all nesting 
levels. Below we first introduce notations and explain 
the existing proof for UDD [BHS] using the language that 
can be naturally generalized for QDD/NUDD. 

UDD. — Let us first consider the UDD sequence 0^ 
for a qubit-environment interaction 



H (r) = So®Bo (r) + Si(g)Bi (r) 



(1) 



where Sq — I and — az, and the time-dependent bath 
operators are analytic with series expansion Ba (t) = 
J2'^=o^o',p^^ for a = 0,1. The UDD sequence uses N 
TT-pulses (i.e., ax rotations) applied at times t\ = TA\, 
where 



Aa = sin^ 



2{N + 1) 



(2) 



for A = 1, 2, • • • ,N. (An extra pulse is required at 
tn+1 ^ T tor N odd [Hj.) It is convenient to consider 
the toggling frame associated with the qubit. In this 
frame the qubit-environment Hamiltonian is modulated 
in time as: 

H (r) - Fo (t/T) So® Bo (t)+Fi (t/T) Si^B^ (r) , (3) 

where (t) = (-1)"'^ for t G (Aa, Aa+i]. The unitary 
evolution operator of the toggling frame Hamiltonian is 

U (T) = Texp — i H (t) dr , where T is the time- 
ordering operator. U (T) has Dyson expansion 

CO s 

where <S'(0q,^.) = Sai ■ ■ ■ S^, and the coefficient J-a\'^.'.'. 
can be obtained by the following integral [6] 



dt. 



dh\{F^^{t,)f^\ (4) 



When ©J^ittj = 0, the operator <5'(0qj.) = / is the iden- 
tity operator that acts trivially on the qubit. Hence 
we only need to consider the terms with ©*^iQ;j 7^ 
that act non-trivially on the qubit. To show the uni- 
versality of the UDD sequence, we just need to prove 
U {T) ^ I ®Ub (T) -I- O (r^+i), which can be reduced 
to verifying 







(5) 



for ©l^iCtj ^ and s + X]^=iPj < ^- Eq-([5|) resem- 
bles the proof of Ref. [Slfor universality of UDETT^The key 
difference is that here additional indices {oij} are intro- 
duced to label different possible qubit operators (/ and 
Gz) which will be necessary for the proof of universality 
of QDD and NUDD. 



QDD. — Let us now consider the QDD sequence [T3] 
for generic interaction between a single qubit and envi- 
ronment, 



(6) 



where Sa — IjCrxjCTyTcrz for binary vector labels a — 
(02, ai) = (0,0), (1,0), (1,1), (0,1), respectively.^ For 
Pauli matrices, one can verify that SaSa' — iSa®a', 
where © represents pairwise binary addition without 
carry (e.g., (0,1) © (0,1) = (0,0)). The QDD sequence 
consists of two nesting levels of UDD with a total of 
Q — (A^ + 1)^ pulses [E]. The pulses cTx and (Tz are 
associated with the first and second levels, respectively. 
To label these Q pulses, we introduce the vector la- 
bel A = {I2, h) G {O,--- ,N} ® {!,■■■ ,N + 1} with 
I2 (N + l)+li e {!,■■■ ,Q} m]. The Ath pulse is applied 
at time t\ = TAf^i^ i^-j with 
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h+i 



AiJAz,. 



(7) 



The toggling frame Hamiltonian for the QDD sequence 



(8) 



where F^ {t) = (-1) 
One can verify that 



for t e (A(z^,,^),A(i^j^+i)]. 



F^ {t/T) F^, {t/T) = F^^^. {t/T) , 



(9) 



which will be useful for our proof of universality of the 
QDD sequence. Using the Dyson expansion of the uni- 
tary evolution operator of H (r) for QDD, we obtain that 
to show the suppression of both the dephasing and bit- 
flip errors up to O {T^'^^^ for small T, it is sufficient to 
prove Eq.([5]) for ®]^^aj 7^ (0,0) and s + Y.]=iPj < N. 
This is very similar to UDD, the difference being that aj 
is now a two-component binary vector. 

NUDD. — The NUDD sequence is a generalization 
of QDD from one-qubit to multiqubit systems [TBl- 
I18| . For m-qubit system, the most general system- 
environment interaction can be written as Eq.([6]) 



and a 



with Sa = (JyJ (g) cri^_i (8) • • • (X) (Til 
{0,2 ) e {0, l}®^"" for all generators. The 

Pauli operator of the jth qubit is ay/ = 1, ax ,cry \az 
for (a2j, a2j-i) = (0, 0) , (1, 0) , (1, 1) , (0, 1), respectively. 
The NUDD sequence consists of 2m nesting levels and 
a total of Q — {N + 1)^™ pulses. The decoupling 
pulse is ai/^ for (2j - l)th level and ai^^ for 2jth 
level. Similar to QDD, we introduce the label A = 
l2m^-^,■■■ Ji) e {O,--- ,iV}^""-'(8){l,--- ,N + 1} 

Q}. The Ath pulse is 



{hm, l'2m-l, 

with {N + i)''-' e {1, 

applied at time ta = TAf^i^^ . 
cursively 



which is defined re- 



A 



A,^ + (A,^+i-AiJA(,_^, 



(10) 
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for r = 2, • • • , 2m. As illustrated in Fig.[l] the timing for 
the pulses has a self-similar structure |25j. The toggling 
frame Hamiltonian for the NUDD sequence is the same as 
Eq.Q, with a summed over all 4™ generators. Similar 
to UDD and QDD, to show universality of the NUDD 
sequence, we just need to prove Eq.([5| for ^ 

and s + — ^ ■ We can view UDD and QDD as 

special cases of NUDD with one and two nesting levels, 
respectively. Since the universality of UDD, QDD, and 
NUDD all relies on verifying Eq. ([5| , we will give a general 
proof of Eq.([5| in the rest of the paper. 

Universality Proof. — To prove Eq.([5]), we represent 
each integration over tj as a linear operator acting on 
functions of tj, which generates a function of ij+i. Thus, 
the process of multiple integrations can be thought of 
as a discrete quantum walk in a functional space. We 
choose the basis of this functional space according to the 
following consideration — the functional basis should be 
complete with respect to the operation dt' Fa (f) t^ 
for all relevant orders up to O (r^+^). Since Fa (t) is a 

piece- wise analytic function with at most Q = (iV -I- 1)^™ 
discontinuities, it will be convenient to use piece- wise an- 
alytic functions as our functional basis. In addition, we 
only need to consider all polynomials up to power iV -t- 1 
to characterize all the effects up to O (r^+^) . Therefore, 
we choose basis that consists of (iV + 1) • Q functions 



where «i < -"^-l + Ei^iPj <N-l,/3o 

and c"^.'"^ij are possibly non-zero coefficients. The 

Q-vectors ul) and \uft,) are determined by Jq dts and 1 
in Eq. (12), respectively, with vector elements \ul)-^ = 
A>,+i — and \ufj}^ = 1. The Q x Q matrices are 



{Dp)l = A,<5f - (A,+i - A,) (-1 

A"=A+1 



(13) 

4',. 



Discrete Quantum Walk. — The last step is to verify 
{uL I DX Dl-_\ ■ ■ ■ \Bp,UR} = (14) 

for (3o ^ and ij < N — 1. The left hand side is 

the amplitude of a discrete "quantum walk" from initial 
state \ul) to a target state \B/3gUii) in the functional 
basis. Each multiplication of Dp corresponds to one step 
of a quantum walk. We need to show that the target state 
amplitude is zero when the number of steps is X]j=i *j — 
iV- 1. 

To understand the quantum walk we choose a conve- 
nient basis 



Xk {t) = ^ C„,A??A (t) 



(15) 



(11) 



with r]\ (t) = 1 for i e A^+i] and r]\ {t) = other- 
wise. Here q € {0, 1, • • • , N} is the polynomial label and 
A e {0, • • • , N}'^^"' is the pulse label. 

Then, we use Eq.Q to rewrite the integral 



*^ ai ,as 



dtstP' X 



dt^^it 



s-ll-s-l 



dtit\' X Fp„ {ti) X 1, 



(12) 



where /J dt' = Fp (t) /J dt'Fp {t') and 
For s + X]j=i?'i — '^^^ compute the integral 

Fa\','.-- 'a° using functional basis {r]q \ (t)}. For each oper- 
ation O, let us define the matrix form as O • rjg^x = 

x Vq'.x', where summation over repeating indices is 
implied. For example, the multiplication t ■ r]q^\ — »7g+i,A 
has the matrix form M."^'-^ = S^^iS^^ , with Kronecker 
delta function 6y. The other operations are listed in Ta- 
ble H 

Using block diagonal properties of matrices involved, 
■Fal'..-- can be further reduced as multiplication of Q x 
Q sub-matrices Bp and Dp [21 



VL ■ Mf= • Gp^_ 

il,--- ,i 



MP-i Gp,-MP'- Fp„vl 

{uL\DlDl-_\---DlBpJun), 



with orthogonal transformation 

■ (2/^ + 1) fc^ + (iV - 



Ck.A 



n 



1) {kr-l 



N + 1 



(16) 

where k = (fem, fc2m-i, ■ • • ^^i) e {0, ■ ■ • ,iV-l-l}'^^" 
and fco = 1 [H]. In the basis {\k)} with |k) = Xk, 
the initial state is \ul) oc |(1,--- ,1)). Since /Jg ^ 
(with br* = 1 and &,.<r* = 0), the target state is 
\BpqUji) — X]k# ""(So-fi* 1'^*) where the r*th element of 
k"^ is k% = N + 1. After some calculation [21 , it can be 

shown that Dp \k) = Y^lltl ' ' ' Y.klZ=o ^^^'^^^ I'*') where 
dK',K denotes possibly non-zero coefficients. Since each 
step of quantum walk only increases the index kr by at 
most 1 unit, it requires at least TV steps to walk from 



Operation 


Matrix/ Vector Form 


t ■ VqA = Vq+1,\ 




Ff3{t)-ll,.^ = {~lt^nq,X 




Il^dt' ■ rig.x = {Gfs)l[f riq'.y 


w,x' 

y^f^lq.X - q + 1 


j'odt-rig^x = (vl)^^^ 




1 = ("fl)<,,A'?9,A 


i^n),,,=S^, 



{i,>0} 



TABLE I: Matrix/vector forms of operations. Here Sy is the 
Kronecker delta function. Bp and Dp are Q x Q matrices as 
defined in Eq. ([T3|. 
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FIG. 2: (color online). Evolution under discrete quantum walks for A'^ = 4. The symbols represents different types of states: 'S' 
for the initial state, ' V' for states explored by a quantum walk, for unexplored target states, '@' for explored target states, 
'.' for remaining unexplored functional basis states. At least A'' = 4 steps are necessary to reach the target states, which is an 
illustration of Eq. (141. (a) For UDD, a discrete quantum walk occurs in a one-dimensional functional basis, with initial state 
target state \N + 1), and quantum walk matrix Di. (b) For QDD, a discrete quantum walk occurs in a two-dimensional 
functional basis, with initial state target states {]«:*)}, and a quantum walk matrix -D(oi). 



kr* = 1 to kf, — N + 1. Therefore, there is zero am- 
phtude in the target states when the number of steps 



J2j=i h — N—l. This completes the proof of Eq.( 14 1 and 
imphes the universality of NUDD. We emphasize xKat af- 
ter the basis in Eq. (11 1 is introduced, all operations are 



matrix multiplications, and all our analytical statements 
have been explicitly checked numerically. 

We can illustrate the quantum walk for special 
cases: (1) For UDD, the orthogonal transformation 

Ck,\ = (^{ki).Xh) = ^iTv+rf simply the Fourier 
transformation as in Ref. 6, The functional basis 
{\k)} with K — (fci) £ {I,-- - ,^+1} forms a one- 
dimensional array. As illustrated in Fig. [2^, it will 
take at least N steps to walk from k = (1) to k'^ — 
(A^ + l). (2) For QDD, the orthogonal transformation 
is a little more complicated c„,a = Hk2,ki),(h,h) = 



(-1) 



(fci-l)/2_ 



-I \ ki /2 

^1 ' cos 



J, 2i2 + l TT 

'^2 N+i 2 

L 2I2 + I TT 

^2 N+1 2 



sm 



I, 2ii+l TT 

'^1 N+1 2 

J. 2;i+i TT 

N+1 2 



for odd ki 
for even ki. 



and (- 

The functional basis with n — (fc2,fci) forms a 

two-dimensional array. As illustrated in Fig. ^p, it will 
also take at least N steps to walk from k — (1,1) to 
K* = {N + l,ki) or (fc2,A^+l). (3) For NUDD, the 
functional basis forms a 2TO-dimensional array. Similar 
to UDD and QDD, it will take at least N steps to walk 
from K — {!,■ ■ ■ , 1) to with kf = N + 1. 

When the suppression order is N,, for the rth nesting 



level of NUDD, the overall suppression of decoherence is 
O (T^ limited by the lowest suppression order N* = 
min[iVi,--- ,N2rn]- 

Summary & Outlook. — We have proved the univer- 
sality of the NUDD sequence, which can restore an un- 
known initial state of m-qubit system to O (^T^+^) us- 
ing 2m nesting levels and {N + l)^™ pulses, indepen- 
dent of the details of the system-environment interac- 
tion. The NUDD sequence is experimentally feasible, 
because it requires only poly (N) pulses acting on indi- 
vidual qubits. Our work illustrates a general connection 
between DD problems and discrete quantum walks. The 
techniques developed can be used to address a variety 
of questions, such as investigation of environment corre- 
lations using DD, schemes of efficient DD for particular 
system-environment interaction, and the combination of 
multiqubit DD schemes with quantum error correcting 
codes I22J and quantum algorithms. 

Note added. — After submission of this work, a 
preprint of closely related work by Kuo and Lidar [53] 
became available, with a different proof of the universal- 
ity of QDD. 
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